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We report on the observation of anti-ferromagnetic correlations of ultracold fermions in a variety 
of optical lattice geometries that are well described by the Hubbard model, including dimers, ID 
chains, ladders, isolated and coupled honeycomb planes, as well as square and cubic lattices. The 
dependence of the strength of spin correlations on the specific geometry is experimentally studied 
by measuring the correlations along different lattice tunneling links, where a redistribution of corre¬ 
lations between the different lattice links is observed. By measuring the correlations in a crossover 
between distinct geometries, we demonstrate an effective reduction of the dimensionality for our 
atom numbers and temperatures. We also investigate the formation and redistribution time of spin 
correlations by dynamically changing the lattice geometry and studying the time-evolution of the 
system. Timescales ranging from a sudden quench of the lattice geometry to an adiabatic evolution 
are probed. 
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Understanding the mechanisms underlying quantum 
magnetism is among the most thought-provoking chal¬ 
lenges of quantum many-body physics. At the center of 
these efforts is the interplay between the emergence of 
magnetic correlations and the underlying lattice geom¬ 
etry [1]. Extensive research has been carried out using 
different materials, as well as theoretical and numerical 
methods, which enhanced the understanding and also 
triggered unforeseen questions [2, 3]. A new develop¬ 
ment is the study of quantum magnetism using ultracold 
fermionic atoms in optical lattices. This technique of¬ 
fers a uniquely direct link between experimental obser¬ 
vations and theoretical models, a key element for quan¬ 
tum simulation [4, 5]. It also promises unprecedented 
dynamic control over lattice parameters and geometry [6- 
10], which can give an entirely new perspective on out-of¬ 
equilibrium properties of quantum spin systems [11]. In¬ 
deed, anti-ferromagnetic spin correlations were recently 
observed in the Hubbard regime of an optical lattice, first 
in an anisotropic [12] and later in an isotropic simple 
cubic geometry [13]. In both situations highly sophisti¬ 
cated computational methods were required for compar¬ 
ison between experiment and theory [13-16]. Along a 
different line, making use of ultracold bosons, progress 
was made in simulating static and dynamic properties of 
classical and quantum mechanical spin models in theo¬ 
retically more tractable regimes [17-23]. 

In this letter we explore the emergence of anti¬ 
ferromagnetic spin correlations in different lattice geome¬ 
tries of varying dimensionality, also including crossover 
configurations between different geometries. The dy¬ 
namic control over the geometries enables us to study 
the formation dynamics and the redistribution time of 
spin-correlations, where the explored timescales range 
from the sudden to the adiabatic regime. The starting 
point of the experiment is a harmonically confined ul¬ 


tracold Eermi gas of in a balanced two-component 
spin-mixture with repulsive interactions [12]. The atoms 
are prepared in the two magnetic sub-levels mp = —9/2, 
—7/2 of the E = 9/2 hyperfine manifold, and the s-wave 
scattering length is tuned between 136.4(5) — 149.0(3) ao 
via the Eeshbach resonance located at 202.1 G (ao de¬ 
notes the Bohr radius). Eor all experiments, the atom 
number is 140(30) x 10^ with 10% systematic error and 
the temperature is 0.09(1) Tp, where Tp is the Eermi tem¬ 
perature. After the preparation, the atoms are loaded 
into the lowest band of a tunable-geometry optical lat¬ 
tice using an S-shaped ramp lasting 100 ms. The lat¬ 
tice consists of several retro-reflected interfering and non¬ 
interfering laser beams of wavelength A = 1064 nm, which 
gives access to a broad variety of lattice geometries [8, 24] . 
Additionally, in all measurements a 3D harmonic confine¬ 
ment is present in the experiment with a geometric mean 
trapping frequency of u/27r = 57(1) Hz. 

Our experiments are well described by the Eermi- 
Hubbard Hamiltonian 

= - E + h.c.) + E + E ’ 

{ij),a i i,cr 

( 1 ) 

with tunneling tij between nearest neighbours (ij) and 
repulsive on-site interaction U. Here and Qct denote 
the fermionic creation and annihilation operators for the 
two spin states a G Ihe density operator is de¬ 

noted by hia- = cj^Cia- and Vi is the trap energy. The dif¬ 
ferent lattice geometries are realized experimentally by 
independently adjusting the specific values of tij for each 
of the six nearest neighbour links per lattice site of an 
underlying simple cubic lattice. Their strength is con¬ 
trolled via the power of the lattice laser beams [24] . In all 
measurements presented in the following, we adjust the 
scattering length with a Eeshbach resonance such that 
U/h = 0.87(2) kHz and set the total bandwidth for non- 
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Figure 1. Experimental observation of the dependence of anti¬ 
ferromagnetic correlations on lattice geometry. The trap av¬ 
eraged correlator —{SfSf_^i) — is measured along 

the strong links for various lattice geometries, which differ in 
the number of strong nearest neighbour links Z. A schematic 
view of the lattice geometries is additionally shown, where 
the strong tunneling links are indicated by a bar. For all data 
points the bandwidth Wjh — 2.6(1) kHz and the on-site in¬ 
teraction Ujh — 0.87(2) kHz is constant. Error bars denote 
the standard error of 50 measurements. 


interacting particles to Wjh = 2.6(1) kHz [25]. 

After loading the atoms into the desired lattice ge¬ 
ometry, we measure the trap-averaged magnetic spin 
correlations emerging on neighbouring sites in the low- 
temperature many-body state of the quantum gas. Our 
detection is similar to the method used in previous ex¬ 
periments [12], and is presented in full detail in the Sup¬ 
plemental Information [24]. The spin correlations are 
measured on every second lattice link, between nearest 
neighbours % and i +1, and along the transverse spin axis 

0,i+i = (2) 

Here denote the standard spin vector operators 

for a spin-1/2 system on site i, and (...) denotes the 
trap average. For SU(2) symmetry, is equal to 

The detection protocol allows us to mea¬ 
sure both anti-ferromagnetic and ferromagnetic configu¬ 
rations, corresponding to positive and negative values of 
respectively. 

In a first measurement, we investigate the strength of 


spin correlations in several different lattice geometries. 
Starting from an underlying simple cubic lattice, the tun¬ 
neling is enhanced along Z nearest neighbour links and 
takes the value whereas the tunneling along the re¬ 
maining 6 — Z links is ts/5. The geometries realized in 
this manner are, sorted by increasing number of strong 
nearest neighbour links Z: dimerized, ID chains, honey¬ 
comb planes, ladders, square, coupled honeycomb planes, 
and cubic (see Fig. 1). We measure the correlations 
along the strong tunneling links. As shown in Fig. 1, the 
strength of the correlations depends on the specific geom¬ 
etry with values ranging between 0.084(1) and 0.010(1) 
for the trap averaged value, and is generally smaller for a 
larger number of strong tunneling links. In the isotropic 
cubic lattice as well, which has the largest value of Z 
(Z=6), we detect anti-ferromagnetic correlations in the 
system [13]. 

The observed dependence of the spin correlator on Z 
can be understood with two simple arguments in a ho¬ 
mogeneous system. Owing to the isolated nature of the 
system, the total entropy rather than the temperature is 
constant for different Z. First, with a finite entropy, the 
presence of two different energy scales associated with 
different tunnelings directly affects the magnetic corre¬ 
lations and leads to a redistribution of spin correlations 
between the strong and the weak links. For a large num¬ 
ber of weaker tunneling links, more low-energy states are 
accessible. Thus, a finite entropy mainly leads to thermal 
fluctuations within theses states, and the magnetic corre¬ 
lator on the strong links is high. However, if the number 
of weak links decreases, thermal fluctuations along the 
weak links alone are not sufficient to account for the total 
entropy, and additional thermal fluctuations are also dis¬ 
tributed on the strong links, therefore reducing the strong 
link correlator. Second, even at zero entropy, quantum 
fluctuations play a significant role. In lower dimensions 
this generally leads to enhanced short-range spin corre¬ 
lations [27]. In both cases, the correlator is expected to 
decrease as Z is increased. This is in accordance with 
previous measurements in the specific cases of a dimer¬ 
ized lattices [Z — 1) and ID chains [Z = 2) [12, 15]. 

While these two effects predict a dependence only on 
Z, the lattice geometry itself (for the same value of Z) 
will also affect the strength of the spin correlations, most 
importantly at low temperatures. In fact, in the limit 
of vanishing temperatures the state of the system and 
its phase diagram will be entirely determined by the in¬ 
terplay between geometry and magnetic ordering. The 
different values observed for the two Z = 3 geometries 
(ladder and honeycomb planes) suggest that effects of 
the lattice geometry are already starting to play a role at 
the temperatures reached in the experiment. An analysis 
of the exact extent of this effect may be achieved by a 
detailed experiment-theory comparison, which is beyond 
the scope of this work. 

In addition to the aforementioned contributions, the 
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Figure 2. Spin correlations for a crossover between two differ¬ 
ent lattice geometries, (a) Scan between a square (Z = 4) and 
a ID chain (Z = 2) geometry and (b) scan between a ID chain 
(Z = 2) and a dimerized (Z = 1) geometry. In both cases, the 
strong to weak tunneling ratio ts/t varies between 1 and 5. 
Blue and red symbols denote the measured spin correlations 
along the strong and weak links respectively. Gray symbols 
show their sum. The error bar on the tunneling ratio denotes 
the uncertainty on the lattice parameters, while the data is 
the mean di the standard error of at least 25 measurements. 


underlying harmonic confinement also plays a central role 
for the value of the trap-averaged correlator. Within the 
local density approximation, both the chemical potential 
at the center of the trap and the temperature are deter¬ 
mined by the total atom number and entropy. As the 
equation of state of the system depends on the lattice 
geometry and Z, both the density and entropy distribu¬ 
tion change with the geometry, which directly affects the 
magnetic correlator. 

To further study the impact of geometry on magnetic 
correlations, we measure their strength for a crossover 
regime between two lattice geometries with different 
numbers of strong links Z. In the experiment geome¬ 
tries with a different Z can be smoothly connected by 
adjusting the strength of the individual tunneling links 
[24]. We scan between a square (Z = 4) and ID chain 
(Z = 2) geometry by choosing the strong-to-weak tun¬ 
neling ratio ts/t in the range 1 — 5, while keeping the 
tunneling in the orthogonal direction at see Fig. 

2a. A second scan between a ID chain and dimerized 
(Z = 1) geometry is shown in Fig. 2b. We measure spin 
correlations either along the strong or the weak tunnel¬ 
ing links. In both cases, correlations on each link start 
from the same value. As ts/t is increased, the correla¬ 


tions along the strong links are enhanced, whereas the 
correlations along the weak links decrease. Interestingly, 
the correlations change more rapidly with increasing ts/t 
for the scan in Fig. 2a as compared to Fig. 2b, which is a 
consequence of the underlying lattice geometry. For the 
final configuration ts/t = 5, the correlations on the weak 
link have nearly vanished, whereas the correlations on the 
strong links have saturated at a high value. This indi¬ 
cates that the thermal fluctuations occur predominantly 
on the weak links. Consequently, the weaker couplings 
can be neglected in this case, and the dimensionality of 
the lattice is effectively reduced for our total entropy and 
atom number. Interestingly, the reduction of the dimen¬ 
sionality occurs at different ratios of ts/t, depending on 
which geometry is considered. 

These measurements demonstrate that spin correla¬ 
tions redistribute between the strong and weak links 
when changing Z. Yet, this does not necessarily imply 
that the sum of spin correlations is constant. We find the 
sum of correlations to be approximately constant in the 
scan of Fig. 2a, whereas it increases significantly with 
dimerization in the scan of Fig. 2b. This observation 
might be related to the opening of a finite energy gap 
in the energy spectrum for a strongly dimerized lattice, 
which causes entropy redistribution within the trapped 
system and enhances the overall spin correlation strength 
[ 12 ]. 

The tunability of our lattice also allows us to experi¬ 
mentally measure the timescales for the formation and re¬ 
distribution of spin correlations when dynamically chang¬ 
ing the lattice geometry. This is done by measuring the 
strength of spin correlations when the lattice geometry is 
changed on variable timescales. For simplicity we start 
with a ramp where the initial and final lattice geometry 
are the same up to a rotation: starting from a ID chain 
geometry we ramp via a square lattice to a ID chain 
lattice again, but with strong tunneling along the per¬ 
pendicular direction. We always include an additional 
wait time before the ramp such that the total time in the 
optical lattice is constant, see Fig. 3a. The spin correla¬ 
tions are measured immediately after the ramp along the 
two different directions. 

The observed dependence of the spin correlations on 
the total ramp time r is shown in Fig. 3a. For fast 
ramps, r < 1ms, the spin correlations remain unchanged, 
at the values without ramping. Here, a non-equilibrium 
state is formed with several charge and spin excitations, 
which decay when allowing for an additional wait time 
after the ramp - thus changing the detected value of spin 
correlations (open symbols in plot Fig. 3a). On interme¬ 
diate ramp times, r ^ 4 ms, the spin correlations change 
symmetrically along the two directions. This timescale 
is comparable to the underlying tunneling time between 
h/ts = 2 ms and h/t = 10 ms during the ramp. For very 
slow ramp times we observe within error bars a 100% 
transfer of spin correlations from the previously strong 
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Figure 3. Dynamics of spin correlations. The lattice is ramped within a time r a from a ID chain geometry to the same ID 
chain geometry rotated by 90° and b from a dimerized to ID chain geometry [24]. Blue and red points denote the measured 
spin correlations along the previously and new strong links. For all closed symbols the correlations are measured immediately 
after the ramp, whereas open symbols include an additional wait time of 50 ms. The gray solid lines indicate the reference 
value for the correlations in the initial geometry without lattice ramp, and the shading denotes the error on this value. Error 
bars are the standard error of at least 25 measurements. 


to the new strong links. When waiting 100 ms in this 
case, the magnetic correlations decrease (most likely due 
to underlying heating of the gas), but agree with the 
value when loading directly from the harmonic trap to 
the final lattice geometry and waiting for the same total 
time. These observations are in agreement with a fully 
adiabatic ramp to an equilibrated final state for the slow¬ 
est ramp times, as the initial and final lattice geometries 
are the same with the same density and entropy distri¬ 
butions. 

The situation changes considerably for a ramp with 
different initial and final lattice geometries. Here, we 
start from a dimerized lattice with a ratio of tg/t = 5 
for adjacent tunneling links and ramp to a ID geome¬ 
try without dimerization. Immediately after the lattice 
ramp we measure the spin correlations on the initially 
strong and weak links along the ID direction, see Fig. 
3b. As before, for fast ramps the spin correlations can¬ 
not redistribute and are nearly unchanged as compared 
to the case without ramp. When adding an additional 
wait time after the fastest ramp, the spin correlations 
change, signaling a decay of the created excitations in 
this case. The behavior is different for intermediate ramp 
times: while the correlations on the initially strong links 
decrease very quickly, slower ramp times r ^ 1 ms are 
necessary for the correlations on the initially weak links 
to change. This may originate from the difference in the 
overall tunneling timescale during the ramp for the two 
links. For the slowest ramps, the correlators along the 


original strong and weak link are identical. With an ad¬ 
ditional wait time, they both decay to the same value, 
again owing to heating. This indicates a final state close 
to equilibrium. In contrast to the previous measurement, 
the gap between the ground and excited states closes dur¬ 
ing the ramp, since the singlet-triplet gap has vanished 
in the non-dimerized geometry. Yet, the observed spin 
correlation value agrees with a reference when loading 
directly into the final lattice geometry and holding the 
remaining time. Consequently, ramp times correspond¬ 
ing to a few tunneling times are already sufficient to reach 
equilibrium. 

In conclusion, we have observed anti-ferromagnetic 
correlations in a variety of lattice geometries and stud¬ 
ied the redistribution of correlations between strong and 
weak links. Extending our work to lower tempera¬ 
tures will allow addressing open questions on the low- 
temperature phase diagram of the Hubbard model in 
complex lattice geometries. In combination with re¬ 
cent efforts on creating spin-dependent optical lattices 
with very low underlying heating [28], the anisotropic 
XXZ Heisenberg model and the Falikov-Kimball model 
could be realized and studied experimentally with ul¬ 
tracold fermions [29]. Our results on quantum spin dy¬ 
namics demonstrate that ultracold fermions in optical 
lattices are well suited to study open questions in out- 
of-equilibrium many-body spin systems, where theoret¬ 
ical methods become extraordinarily difficult [30, 31]. 
Furthermore, the observed rapid formation timescales of 
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spin correlations offer very promising perspectives for the 
implementation of sophisticated entropy redistribution 
schemes based on trap shaping and dynamically changing 
lattice geometries, which are expected to result in overall 
lower temperatures [32-34]. 
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SUPPLEMENTAL MATERIAL 
Preparation and optical lattice 

For our measurements we prepare a quantum degen¬ 
erate cloud of K fermions, harmonically confined in 
a crossed optical dipole trap. This is done by evapora- 
tively cooling a balanced spin mixture of the two mag¬ 
netic sub-levels mp = —9/2, —7/2 of the F = 9/2 hy- 
perfine manifold. We reach temperatures of 0.09(1) Tp, 
where Tp is the Fermi temperature. This corresponds 
to an entropy of 0.87/^^ per particle in the harmonic 
trap. The atom number is set to a value of 140(30) x 10^ 
for all measurements. The repulsive on-site interaction 
energy U is tuned via a Feshbach resonance located at 
202.1 G by changing the s-wave scattering length between 
a = [136.4(5) — 149.0(3)] ao, where ao denotes the Bohr 
radius. 

After the preparation, the fermions are loaded into a 
tunable-geometry optical lattice within 100 ms using an 
S-shaped ramp. The lattice potential is created by sev¬ 
eral retro-refiected interfering and non-interfering laser 
beams of wavelength A = 1064 nm, which gives access 
to a broad variety of lattice geometries described by the 
following potential [SI]: 

V{x, y, z) = —V^cos^{kx + 0/2) — Vx cos‘^{kx) 

— UyCOS^(%) — cos‘^ {kz) 

—2a\/VxVY cos{kx) cos{ky) cos (p^ (SI) 

where xyz lattice depths of each single beam 

in direction x^y^z given in units of the recoil energy 
Er = h^/2mA^, k = 27r/A, visibility a = 0.81(1) and h 
denotes the Planck constant. We set 0 = ir and interfer- 
ometrically stabilize the phase ip = 0.00(3)7r. The lattice 
depths XYZ independently calibrated using Ra- 
manNath diffraction on a ^^Rb BoseEinstein condensate. 

Each lattice geometry shown in Eig. 1 is obtained by 
individually adjusting the lattice depths Vx xyz' 
lowing this procedure the connectivity Z (i.e. the num¬ 
ber of strong nearest neighbor links) can be varied from 1 
(dimerized lattice) to 6 (cubic lattice). This is achieved 
by enhancing the tunneling rate of Z nearest-neighbor 
lattice bonds, with respect to the remaining 6 — Z bonds 
per site. Due to the tunability of our optical lattice we 
can realize the seven individual lattice geometries shown 
in Eig. 1, as well as the crossover configurations, such as 
the ones shown in Eig. 2. Additionally it allows us to 
ramp from one lattice geometry to another on a variably 
time-scale, which is used for the dynamics measurements 
shown in Eig. 3. 

To set the correct lattice depths we calculate all Hub¬ 
bard parameters numerically from the lattice potential 
using Wannier functions. Our computational method is 
based on their definition as eigenstates of band-projected 
position operators [S2]. Erom the overlap integrals of the 


Wannier functions we obtain the tunneling tij between 
nearest neighbors and the interaction energy U. Eor all 
our measurements, we choose the lattice depths in such a 
way that the total bandwidth of the single orbital tight- 
binding model is always W/h = 2.60(6) kHz, independent 
of the particular lattice geometry. The expression is given 
by VF = 2 toi, where i denotes the sum over all 6 near¬ 
est neighbors of any given lattice site in the underlying 
simple cubic lattice. Since the Wannier functions depend 
on the exact geometry and lattice depths, we have to ad¬ 
just the scattering length by tuning the magnetic field 
of the Eeshbach resonance to ensure a constant on-site 
interaction U/h = 867(15) Hz. The exact parameters for 
all the lattice geometries with different connectivities re¬ 
alized in the experiment are summarized in Table SI. 

Due to the changing harmonic confinement as a result 
of the intensity variation of the lattice beams, we have 
to adjust the optical power of the additional dipole trap 
to set a constant geometric mean trapping frequency of 
uj/2'K = 57(1) Hz for all measurements. However, the 
individual trapping frequencies in x-, y- and 2 :-direction 
of the harmonic confinement are changing for each lattice 
geometry. 

Detection scheme for spin-correlations 

The detection scheme for nearest-neighbor anti¬ 
ferromagnetic spin correlations consists of several steps 
(see Eig. SI) and closely follows the procedure used in 
previous experiments [S3]. After preparing the many- 
body state in the chosen lattice geometry, we perform 
a sudden ramp of the lattice depth to freeze out the 
atomic motion. Eor this we always ramp into a deep 
simple cubic lattice with V^^y^ = [30,0,30,35] 
within 100 /is independent of the starting geometry of 
the optical lattice (i.e. Z=1..6). The overall mea¬ 

surement projects the state of the initial lattice onto 
the singlet \s) = (It,!) — \iA))/V^ and triplet states 
l^o) = (It^i) + \iA))/^ oa Ihe respective set of links 
where the spin correlations are measured (see Eig. Sib). 
Erom this we determine the fraction of atoms forming 
singlets Ps and triplets pto oa neighboring sites. An im¬ 
balance in these two fractions corresponds to magnetic 
ordering on neighboring sites. The population imbalance 
is proportional to the transverse spin correlator [S3] 

Q,i+1 = -{SfY+i) - = iPs - Pu)/2. (S2) 

which is measured on every second lattice link with near¬ 
est neighbors i and i + 1 and along the transverse spin 
axis (see Eig. Sib). Due to our detection process all 
quantities are averaged over the entire atomic cloud. As 
doubly occupied sites do not contribute to the spin corre¬ 
lator and hinder the detection process, they are removed 
in the deep cubic lattice using a series of Landau-Zener 
transfers before detecting the singlet and triplet fractions. 
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Figure SI. Probing nearest neighbor spin correlations, (a) The detection scheme of the anti-ferromagnetic correlations contains 
several steps and follows closely the procedure used in previous experiments [S3] . After preparing the system in a thermalized 
many-body state with varying geometry, we freeze and convert it into a deep simple cubic lattice. In a second step we apply 
a magnetic field gradient which causes coherent oscillations between the singlet and the triplet state on neighboring sites. In 
a final step we remove the magnetic field gradient and merge pairs of neighboring sites. Depending on the oscillation time we 
detect either the number of singlets or the number of triplets corresponding to the initial many-body state, (b) An overview of 
the set of links between nearest neighbors, along which the spin correlations are measured for the various geometries in Fig. 1 
of the main manuscript, are marked in colored ellipsoids (note that the ladder geometry has been rotated for clarity). 


We measure the singlet and triplet fractions by subse¬ 
quently applying a magnetic field gradient, which lifts 
the energy degeneracy for atoms with opposite spins 
on neighboring sites. This induces a coherent oscilla¬ 
tion between the singlet and the triplet states. The 
singlet-triplet oscillations (STO) have a frequency of 
u = Asto /^5 where Asto is the energy splitting, and are 
only visible if the initial amount of singlets and triplets 
is different, as their respective time evolution is exactly 
out of phase. After a variable oscillation time we remove 
the magnetic field gradient and merge pairs of adjacent 
sites into a single site. We use a 10 ms linear ramp from 


the deep cubic lattice into a deep checkerboard lattice 
{VxxYZ ~ [O 5 25, 30, 35] where tunneling is still 
suppressed and the number of sites is divided by two. 
During merging, the \s) and |to) states are mapped onto 
different bands due to their distinct symmetry of the two- 
particle wave function. The (spatially symmetric) singlet 
state is mapped to two atoms in the lowest band of the 
final state, while the (spatially anti-symmetric) triplet 
state evolves into a final state with one atom in the first 
excited band and one atom in the lowest band. By ad¬ 
justing the oscillation time to a maximum or minimum 
of the STO we can then detect the fraction of atoms in 








































































































the singlet and triplet state by measuring the fraction of 
atoms on doubly occupied sites with both atoms in the 
lowest band. This is accomplished using interaction de¬ 
pendent rf-spectroscopy, which transfers the tuf — —1 
atoms on doubly occupied sites in the lowest band to the 
initially unpopulated ttif = —5/2 state [S4]. 

After all these steps we measure the atomic popula¬ 
tions in the different spin states. We ramp down the 
optical lattice potential and dipole trap within 20ms and 
apply a magnetic field gradient to separate the differ¬ 
ent rriF states in a Stern-Gerlach measurement during 
ballistic expansion. Finally we take an absorption im¬ 
age and apply Gaussian fits to the density distribution of 
each spin component to determine the number of atoms 
in each spin state (exemplary image shown in Fig. SI). 
For the double-occupancy measurement we take into ac¬ 
count the independently calibrated detection efficiency of 
91(3)% and the background offset in double occupancy 
of 2.2% is subtracted (the latter however does not change 
the measured value of the spin correlator). In addition, 
we take into account the separately calibrated decay of 
the singlet-triplet oscillations during the entire detection 
sequence [S3]. 

Calibration of singlet-triplet oscillations 

To test and calibrate the performance of the detec¬ 
tion scheme for measuring spin correlations using singlet- 
triplet oscillations, we prepare an initial state with a large 
number of singlets and very small number of triplets. We 
load an attractive spin-mixture at a = —770(50) ao into a 
deep checkerboard lattice with xyz ~ P ’ 
which leads to a large number of doubly occupied sites. 
We first freeze this geometry and then ramp to repulsive 
interactions with a = 106.5(9) ao- In a second step we 
split the sites of the checkerboard lattice by ramping to 
the same deep cubic lattice which is used in the first step 
of the detection scheme. This converts the double occu¬ 
pied sites to singlets. Although this initial state may be 
far from equilibrium owing to its preparation, it is well 
suited for calibration purposes as it contains a very large 
imbalance in the singlet and triplet fractions, and hence a 
large value for the nearest-neighbor spin correlator. After 
the initial preparation, we apply the entire detection pro¬ 
tocol and measure the double occupancy after merging as 
a function of oscillation time with magnetic field gradi¬ 
ent, see Fig. S2. The fraction of atoms forming singlets 
(triplets) is then given by the maximum (minimum) value 
of the measured double occupancy of the STO, which is 
obtained from a sinusoidal fit to the data. We use this 
measurement to calibrate the magnetic field gradient and 
determine the singlet-triplet oscillation frequency. The 
actual phase of the STO, which originates from residual 
magnetic gradients, is determined by preparing a ther¬ 
mal many-body state in the dimerized geometry of the 



STO time 

Figure S2. Full singlet-triplet oscillation. We measure the 
double occupancy after merging as a function of oscillation 
time during a singlet-triplet oscillation, which is induced by a 
magnetic field gradient in the deep cubic lattice. For this cali¬ 
bration a far-from-equilibrium initial state was prepared with 
a large number of singlets. The observed phase shift is par¬ 
tially caused by a residual slow STO (with a period >100 ms) 
occurring in the deep cubic lattice already before applying the 
magnetic field gradient (originating from residual magnetic 
gradients). Error bars show the standard deviation. 

measurement sequence. From this we determine the max¬ 
imum and minimum oscillation time, which is used for all 
measurements in the main manuscript. 

Fig 1: Connectivity scan 

Using the tunable geometry optical lattice we can indi¬ 
vidually address and adjust the specific tunneling values 
tij for all six links connecting one lattice site. To mea¬ 
sure the dependence of the correlations on the geometry 
we enhance the tunneling G on Z nearest neighbor links 
of the underlying simple cubic lattice with Z ranging be¬ 
tween 1 and 6, while the tunneling on the remaining links 
is set to a much smaller value of ts/5. If all tunneling links 
are set to the same value then ts/h = 217(10) Hz and we 
realize a simple cubic lattice. Table SI summarizes the 
parameters of each of the seven lattice geometries real¬ 
ized for the measurements in Fig. 1. The fraction of 
singlets and triplets is measured on a set of strong links 
for each data point, see Fig. Sib. In the experiment (see 
Fig. 1) we choose for each geometry to always merge the 
same strong links. For the ID chain, square and cubic 
lattice we independently verified that the signal strength 
of the spin correlator does not depend on the exact choice 
of the merged strong link. 

Fig. 2: Crossover between geometries 

In the second set of measurements we explore the 
crossover between two geometries, either from a square 
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Lattice 

Z 

^X,X,Y,Z [^r] 

a [ao] 

U/h [Hz] 

t^/h [Hz] 

t'Jh [Hz] 

ty/h [Hz] 

u/h[nz] 

W/h [kHz] 

Dimer 

1 

[4.925, 0.073, 7.9, 8.19] 

136.4(5) 

866(15) 

635(30) 

136(4) 

130(8) 

130(7) 

2.60(8) 

ID 

2 

[3.205, 0, 9.6, 9.6] 

139.1(5) 

867(15) 

464(12) 

464(12) 

93(6) 

93(6) 

2.60(4) 

2DHC 

3 

[7.545, 0.258, 3.93, 10.7] 

137.5(5) 

868(15) 

352(24) 

75(3) 

363(12) 

72(5) 

2.60(8) 

Ladder 

3 

[7.345, 0.081, 10.45, 4.16] 

136.9(5) 

867(15) 

354(23) 

73(3) 

73(5) 

361(12) 

2.60(6) 

2D 

4 

[4.93, 0, 4.93, 11.6] 

140.9(4) 

867(15) 

295(11) 

295(11) 

295(11) 

59(5) 

2.60(6) 

3DHC 

5 

[9.025, 0.219, 5.37, 5.57] 

143.5(4) 

867(16) 

248(19) 

51(3) 

250(11) 

250(11) 

2.60(8) 

3D 

6 

[6.128, 0, 6.128, 6.128] 

149.0(3) 

867(16) 

217(10) 

217(10) 

217(10) 

217(10) 

2.60(6) 


Table SI. Lattice parameters for each geometry. All lattice geometries of the measurements in Fig. 1 are shown with their 
individual lattice parameters: dimer, ID chains, two-dimensional honeycomb planes (2DHC), two-leg ladder, square (2D), 
coupled honeycomb planes (3DHC) and simple cubic lattice (3D). The different lattice geometries are listed in increasing order 
of the connectivity Z. The graph in Fig. S3 indicates the labeling of tunneling links used in the experiment. To set a constant 
value of the on-site interaction U we vary the scattering length a. With increasing Z the number of strong nearest-neighbor 
links is increasing accordingly (i. e. for the dimerized lattice the strong tunneling is ts/h = 635 Hz, while all other links have 
a much weaker tunneling rate of 130Hz). Due to the exact choice of the individual tunneling links we ensure that the total 
bandwidth W is constant. The respective errors are a result of the systematic uncertainty of our system parameters 



Figure S3. Overview of the different tunneling links in the 
tunable-geometry optical lattice, which are all independently 
tunable. 


lattice {Z = 4) to ID chains (Z = 2), or we start with ID 
chains which are then continuously deformed to a dimer¬ 
ized lattice {Z = 1). In both cases we vary the connec¬ 
tivity by weakening half of the initially strong links to a 
value of t while increasing the other links to tg. At the 
same time the tunneling of the remaining weak links in 
the initial lattice geometry is modified accordingly to 
keep tg/tw = 5 constant. To explain this in more detail 
we consider the case of the ID to dimer crossover. Here 
tx corresponds to the strong tunneling link tg and to t. 
During the crossover as one of the two initially strong 
links in the ID chain (tx and t^) is increased, whereas 
is decreased to a lower value. Thereby the ratio be¬ 
tween tx and can be tuned from 1 to 5. Subsequently 
the tunneling links ty and (corresponding to t^j) are 
increased to keep their ratio with constant. This pro¬ 
cedure allows varying the strong-to-weak tunneling ratio 
tg/t from 1 to 5, thus allowing for a continuous scan be¬ 
tween two geometries. As for the previous measurements 
we also keep the total bandwidth, the geometric mean 
trapping frequency and the on-site interaction constant. 

Each data point corresponds to a series of 25 measure¬ 
ments of the spin correlator after directly loading the 


fermions into the final lattice geometry. Due to the ex¬ 
perimental set-up of the tunable optical lattice it is only 
possible to merge adjacent sites along the x-direction of 
the laboratory frame. In order to measure the spin cor¬ 
relations on neighboring sites along the perpendicular di¬ 
rection as shown in Fig. 2a, the atoms are loaded into a 
lattice geometry that is rotated by 90°, which is achieved 
by exchanging the values for the lattice depths and 
Vy. By merging along the x-direction of the laboratory 
we then effectively merge along the perpendicular direc¬ 
tion of the physical system. Due to the rotation of the 
lattice pattern the effective trapping frequencies along 
all directions of the harmonic confinement are different 
for the two lattice geometries. However, the local den¬ 
sity approximation (LDA), which has been shown to be 
accurate in the temperature regime of our experiment 
[S5], entails that only the geometric mean of the trap 
frequency matters. This implies that the measured re¬ 
sult does not depend on the rotation in the x^-plane at 
constant geometric mean trapping frequency. 

To compare the spin correlations on weak and strong 
links in the dimerized lattice (see Fig. 2b) we need to be 
able to control which of the two sites we merge. This is 
achieved by controlling the phase ip between the X and 
Y lattice beam. During the time spent in the deep lattice 
(when Vx = 0), just before the merging ramp, the phase 
(p is shifted by tt. This ensures that we can merge two 
adjacent sites along the previously weak link. The phase 
shift is generated by an active feed-forward on the phase 
lock using the acousto-optical modulator of the lattice 
beam. 

Fig. 3: Dynamics of spin correlations 

In order to observe dynamics in spin correlations we 
first ramp into the starting lattice geometry, which is ei- 
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ther lattice of ID chain in the x-direction (Z = 2, with 
the same configuration as in Fig. 1) or a dimerized lat¬ 
tice {Z = 1). In contrast to the previous measurements 
we then, in a second step, change the lattice geometry on 
a variable time scale r before detecting the spin correla¬ 
tions. 

In the first measurement shown in Fig. 3a we change 
the ID lattice with strong tunneling links tg initially 
along the x-direction to a geometry with ID chains along 
the perpendicular direction. This is achieved by lin¬ 
early ramping the lattice depths within the ramp time 
T from ~ [3.205,0,9.6,9.6] to Vxxyz ~ 

[9.6, 0, 3.205, 9.6] (at exactly half of the ramp time 
the potential is hence in a square lattice geometry with 
equal tunneling strengths in the x^-plane). The spin cor¬ 
relations are measured in the final geometry, where the 
strong tunneling links are now along the perpendicular 
direction. We do not allow for an additional wait time in 
the final lattice geometry, but directly ramp to the deep 
cubic lattice and follow the detection scheme described in 
the previous section. To keep the time the atoms are ex¬ 
posed to the optical lattice constant, we add a wait time 
before the change in lattice geometry, such that the sum 
of the wait time and the lattice ramp time r is 80 ms. 
To measure correlations along new strong links we fol¬ 
low the same procedure as described in the section of the 


geometry scan. 

In the second set of measurements shown in Fig. 3b 
we analyze the dynamics of spin correlations when lin¬ 
early ramping from a dimerized lattice {Vxxyz ~ 
[4.925, 0.073, 7.9, 8.19] into a lattice with ID chains 
along the x-direction ~ [3.205,0,9.6,9.6] Er). 

We follow the same procedure as before with the only 
difference that the total time in the shallow lattice is 
now set to 100 ms to allow for increased ramp times r. 
To merge on different links (initially we have weak and 
strong tunneling links in the dimerized lattice) we follow 
the same procedure as described for the data in Fig. 2b. 
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